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ABSTRACT 

We report the results from archival XMM-Newton and INTEGRAL observations of 
the Supergiant Fast X-ray Transient (SFXT) IGR J18483-0311 in quiescence. The 
18-60 keV hard X-ray behaviour of the source is presented here for the first time, it 
is characterized by a spectral shape (F ^^2.5) similar to that during outburst activity 
and the lowest measured luminosity level is ^10^'' erg s~^. The 0.5-10 keV lumi- 
nosity state, measured by XMM-Newton during the apastron passage, is about one 
order of magnitude lower and it is reasonably fitted by an absorbed black body model 
yielding parameters consistent with previous measurements. In addition, we find ev- 
idence {r^Z-ba significance) of an emission-like feature at ~3.3 keV in the quiescent 
0.5-10 keV source spectrum. The absence of any known or found systematic effects, 
which could artificially introduce the observed feature, give us confidence about its 
non-instrumental nature. We show that its physical explanation in terms of atomic 
emission line appears unlikely and conversely we attempt to ascribe it to an electron 
cyclotron emission line which would imply a neutron star magnetic field of the order 
of ~3xl0^^ G. Importantly, such direct estimation is in very good agreement with 
that independently inferred by us in the framework of accretion from a spherically 
symmetric stellar wind. If firmly confirmed by future longer X-ray observations, this 
would be the first detection ever of a cyclotron feature in the X-ray spectrum of a 
SFXT, with important implications on theoretical models. 
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1 INTRODUCTION 

IGR J18483— 0311 is a hard X-ray transient discovered 
in outburst with INTEGRAL in 2003 (Chernyakova et al. 
2003). Other hard X-ray outbursts were observed showing 
fluxes and durations up to ~120 mCrab (20-100 keV) and 
~2 days, respectively (Sguera et al. 2007). The X-ray emis- 
sion shows two periodicities: a longer one at ~18.55 days 
which is interpreted as orbital (Levine et al. 2006, Sguera 
et al. 2007), and a shorter one at ~21.0526 s which was 
discovered with the X-ray monitor JEM-X (Sguera et al. 
2007). Giunta et al. (2009) performed a study of the qui- 
escent soft X-ray behaviour of the source (0.5-10 keV) and 
confirmed the pulse period through XMM-Newton data, pro- 
viding a measurement of the spin-period derivative which 
is very likely due to light travel time effects in the binary 
system. A complete monitoring of the X-ray emission from 



IGR J18483— 0311 over an entire orbital period has been 
performed in 2009 with Swift/XRT (Romano et al. 2009). 
The optical counterpart has been identified with a B0.5Ia 
supergiant star located at a distance of 3-4 kpc (Rahoui et 
al. 2008). This, together with the transient nature of the 
source, implies its classification as a Supergiant Fast X-ray 
Transient (SFXTs; Sguera et al. 2005, 2006, Negueruela et 
al. 2005, Sidoh 2009). 



2 XMM-ATiJT^rOJV OBSERVATION AND 
RESULTS 

2.1 Data reduction 

We analyzed an archival XMM-Newton EPIC observation of 
IGR J18483-0311 collected on 2006 October 12 for a total 
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exposure of ~19 ks. However, after rejecting time intervals 
affected by high background the net good exposure time re- 
duced to 14.4 ks. XMM-Newton data were reprocessed using 
version 9.0 of the Science Analysis Software (SAS). Known 
hot, or flickering, pixels and electronic noise were rejected us- 
ing the SAS. The ancillary and response matrices were gen- 
erated using the SAS tasks arfgen and rmfgen. Spectra were 
selected from single and double events only (pattern from 
to 4) for the EPIC pn full frame mode while for both MOS 
cameras patterns from to 12 were selected. Source counts 
were extracted from circular regions of 40 arcseconds radius 
centered on the source for both pn and MOS. Background 
counts were obtained from similar regions offset from the 
source position. The backgrounds do not show any evidence 
for flaring activity, so that the entire nominal exposure times 
were considered. The source net count rates (1-10 keV) 
are the following: 0.097±0.003 counts s"^ (pn), 0.035±0.002 
counts s-i (MOSl), and 0.027±0.001 counts s'^ (M0S2). 
The source is too faint for a meaningful spectral analysis 
with the RGS. Also the spectroscopy with the MOS cam- 
eras does not provide any improvement with respect to the 
analysis of the pn spectrum alone, thus in the following we 
will concentrate only on the source EPIC pn results. In order 
to ensure applicability of the statistic, the extracted spec- 
tra were rebinned such that at least 20 counts per bin were 
present and such that the energy resolution was not over- 
sampled by more than a factor of 3. The spectral analysis 
was performed using Xspec version 12.5; all quoted uncer- 
tainties are given at the 90% confidence level for one single 
parameter of interest. 

2.2 EPIC-pn spectral analysis 

We first fit the 0.5-10 keV EPIC-pn spectrum of ICR 
J18483— 0311 with an absorbed power law model whose best 
fit parameters (NH=7.8+i;? X 10^^ cm"^, r=2.4±0.3) are in 
fair agreement with those already reported by Giunta et 
al. (2009). However, we note that such model gives a sta^ 
tistically inadeguate representation of the observed spec- 
trum being the x^=l.AQ for 52 d.o.f. (see also Giunta et 
al. 2009 where x2 = l-3 for 39 d.o.f.). This motivates us to 
investigate alternative spectral models. An absorbed ther- 
mal bremsstrahlung is still a bad statistical description of 
the data (x^=l-34, 52 d.o.f.), on the contrary an absorbed 
thermal black body results in a much better fit (x^=1.17, 
52 d.o.f.) and yields parameters (N/f=3.4lo 5 x 10^^ cm"^, 
kT=1.35±0.08 keV) consistent with previous measurements 
by Romano et al. (2009) from a Swift /XRT monitoring cam- 
paign covering an entire orbital period. The upper and mid- 
dle panels in Fig. 1 show the absorbed black body fit spec- 
trum and the ratio of data to model, respectively. It is worth 
noting that such spectral fit resulted in a radius of the emit- 
ting black body region equal to ~0.14 km, i.e. consistent 
with a small portion of the neutron star surface such as 
its polar cap region. The unabsorbed (observed) 0.5-10 keV 
fiux is 1.24x10"^^ erg cm"^ s"^ (9.1x10"^^ erg cm"^ s"^) 
which translates into a X-ray luminosity of 1.3x10^^ erg 
s~^ by assuming a distance of the optical counterpart of 
3 kpc. If we follow Sguera et al. (2007) and measure the 
source phase from the epoch of the brightest outburst ob- 
served with INTEGRAL at MJD 53844.2 (phase 0), then 
the XMM-Newton observation took place at orbital phase 
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Figure 1. Upper panel: data-to-model (absorbed black body) of 
the 0.5-10 keV EPIC-pn spectrum of IGR J18483-0311. Middle 
panel: the ratio of data to model is shown for the fit with no 
gaussian line included in the model. Lower panel: the ratio of 
data to model is shown again but this time including a gaussian 
line. 

0.52 (i.e. right during the apastron passage) and the above 
X-ray luminosity value represents the lowest quiescent X-ray 
state of the source ever reported in the literature. For the 
sake of clarity, here and in the following we refer to "qui- 
escence" as to the lowest detectable level of X-ray activity 
from the source (i.e. ~10^^ erg s~^) during which it is still 
accreting matter even if outside its bright outbursts activity 
(i.e. ~10^^ erg s"^). 

2.2.1 An emission-like feature at r^3.3 keV? 

As stated in the previous section, the absorbed black body 
model is a reasonable statistical description of the contin- 
uum; however the ratio of data to model clearly show an ex- 
cess around ^^3.3 keV (see Fig. 1, middle panel), suggesting 
the presence of a possible spectral line. To model such resid- 
uals we added a gaussian emission line with energy, width 
and normalisation free to vary. We obtained a significant im- 
provement of the statistical quality of the fit (x^=0.97, 49 
d.o.f) corresponding to a significance level of confidence of 
99.4% according to a simple F-test (^3cr significance). We 
are aware that the F-test is not a good measure of the actual 
significance of such additional spectral line feature (see Pro- 
tassov et al. 2002), however it could give an indication. To 
this aim, we point out that the obtained low F-test prob- 
ability value should make the detection of the line stable 
against mistakes in the calculation of its significance. Fur- 
thermore, we note that there are in the literature a few cases 
of cyclotron features reported with a F-test probability value 
very similar to our one (e.g. Rea et al. 2003, Ibrahim et al. 
2002). Even using an absorbed power law model, as already 
adopted by Giunta et al. (2009), we still have evidence of 
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the Une although at lower statistical significance (96.7%). 
We measured an energy centroid of 3.28l:°o°oi'i keV; the line 
is unresolved and it has an equivalent width of 148lgQ^ eV. 
The luminosity of the line is ~2.1xl0^^ erg s~^, i.e. ~2% of 
the total 0.5 10 kcV luminosity. The robustness of the line 
detection is attested by the energy ccntroid-normalisation 
parameters confidence contours shown in Fig. 2; moreover 
the line intensity is equal to 5.5^J^ ^ x 10~^ photon cm~^ 
for a significance detection of ^3.5a. Other positive excesses 
in the ratio of data to model arc present in the spectrum at 
~2.5 keV and ~5 keV (see Fig. 1, middle panel), but they 
are not significative. 

In order to rule out an artificial nature of the spectral 
emission line due to instrumental effects, we checked that the 
response matrix did not introduce any strong feature around 
3.3 keV (i.e no sharp variation is present in the spectral 
response). The inspection of the background spectrum did 
not reveal any spectral feature at the same energy. We also 
note that no uncalibratcd instrumental features or edges are 
expected or known close to ~3.3 keV. The absence of any 
known or found systematic effects give us confidence about 
the non-instrumental nature of the observed spectral line. 
In the following we discuss its possible physical origin. 

Firstly, we took into account the possibility that the 
observed emission feature could be due to atomic transition 
lines. In fact, if originating in an accretion gas or corona 
or in the circumstellar material, atomic lines could well be 
present in the X-ray spectra of accreting pulsar and they are 
expected to be mainly narrow and in emission. According to 
the atomic database ATOMDB*, we found that two highly 
ionized Argon atomic lines are expected close to ~3.3 keV 
(i.e. Ar XVIII at 3.32 keV and Ar XVII at 3.14 keV). How- 
ever, we consider unlikely a physical explanation in terms 
of Ar atomic lines because of the following reasons: i) it 
would require an unexplainable high overabundance of Ar 
which is not obvious in the neutron star atmosphere or in 
its surrounding ambient, ii) it would not explain why only 
the highly ionized Ar line is observed and no other lines are 
seen from more abundant elements, iii) it would not justify 
the lack of a similar spectral feature in any of the many 
known accreting X-ray pulsar in HMXBs for which much 
higher signal to noise X-ray spectra are available. Therefore, 
we explored an alternative physical explanation in terms of 
cyclotron line. In fact, depending on the physical condition 
of the emitting/absorbing material in the accretion column, 
cyclotron features from electrons or protons could be ob- 
served in the X-ray spectra of accreting magnetized pulsars 
(Heindl et al. 2004). They provide an important tool for a 
direct measurement of the neutron star magnetic field value, 
since the observed fundamental energy is related to the mag- 
netic field value by E°''J'=11-6(B/10^^ G) in the case of elec- 
trons or by E°fp^=0.63(B/10"G) in the case of protons. If 
interpreted as an electron cyclotron emission feature, the 
putative line observed at ~3.3 keV implies a magnetic field 
of ~2.8xl0" G (if the forming region is situated far above 
the pulsar's polar cap) or alternatively of ^3.7x10^^ G (if it 
is situated close to the neutron star surface, i.e. at the base 
of the accretion column, and so affected by a gravitational 
redshift of z=0.3) . For protons cyclotron feature, the mag- 
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Figure 2. Confidence contours level (68%,95%,99%) for the en- 
ergy line center and normalization parameter values. 

netic field should be even higher by a factor of ~1836 (the 
proton to electron mass ratio) implying a value of ~5xl0^* 
G (i.e a magnetar-type magnetic field). 



3 INTEGRAL OBSERVATIONS AND RESULTS 

The quiescent behaviour of IGR J18483-0311 in the soft X- 
ray band (0.5-10 keV) is fairly known and studied (Romano 
ct al. 2009, Giunta et al. 2009) while above 20 keV is totally 
unknown. With the aim of investigating for the first time the 
hard X-ray quiescent state of the source, we collected INTE- 
GRAL data with IBIS (Ubertini et al. 2003). INTEGRAL 
observations are typically divided into short pointings called 
Science Windows (ScWs) each of ~2000 s duration. The data 
reduction was carried out with the release 7.0 of the Offline 
Scientific Analysis (OSA) software. 

Firstly, we searched the entire IBIS/ISGRI public data 
archive for pointings where IGR .118483—0311 was within 
the fully coded field of view of ISGRI (9° x9°). Subsequently, 
since our aim was to measure the quiescent hard X-ray emis- 
sion, we intentionally excluded those individual ScWs during 
which the source was in outburst, i.e. significantly detected 
at >5cr. By doing so, we collected a total of 401 ScWs which 
were used to generate a mosaic significance map in the 18-60 
keV band for a total exposure of ~ 795 ks. IGR J18483-0311 
was weakly detected at ~7c7 level with an average flux (lumi- 
nosity) of 1.5 mCrab (1.3x10^* erg s^^), likely representing 
the lowest quiescent hard X-ray state of the source. The rel- 
ative spectrum (18-60 keV) is equally well fit using a power 
law (r=2.5±0.8, xl=0-8, 15 d.o.f.) or alternatively a black 
body model (kT=7±2 keV, xl = 0.8, 15 d.o.f.): such spec- 
tral shape is similar to that seen during the outburst activity 
(Sgucra et al. 2007). 

A joint spectral fit to the XMM-Newton and 
IBIS/ISGRI quiescent data, although non simultaneous, was 
also attempted. A multiplicative factor for each instrument 
was included in the fit to take into account the uncertainty 
in the cross calibration of the instruments. We used spec- 
tral models usually adopted to describe the X-ray emis- 
sion from accreting pulsars. Firstly we found that the 0.5- 
60 keV broad band spectrum can be adequately described 
(x^=0.95, 64 d.o.f.) by a phenomenological model such as 
absorbed power law together with a black body at low en- 
ergies, to describe the continuum, plus a gaussian line in 
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Figure 3. Broad-band spectrum (0.5—60 keV) of IGR 
J18483— 0311 fitted with a Bulk Motion Comptonization model 
modified at low energies by photoelectric absorption and a gaus- 
sian line. The lower panel shows the residuals from the fit. 



emission at 3.3 keV (wc loft all parameters free to vary ex- 
cept for the absorption N^f which was fixed to its nominal 
value previously found from the XMM-Newton spectral anal- 
ysis). The best fit parameters are r=2.0tg;^, kT=1.40tg;2^ 
and Ec=3.28±0.05. Subsequently, we used a more physical 
thermal Comptonization model, namely the Bulk Motion 
Comptonization (Titarchuk et al. 1997). The spectral pa- 
rameters of such model are the black body colour tempera- 
ture kTj,6, the spectral index a which is related to the Comp- 
tonization efficiency (when a is smaller then the efficiency 
of energy transfer from the hot electrons to the soft seed 
photons is higher) and finally logA (where A is the illumi- 
nating factor) which gives an indication of the fraction of the 
up-scattered black body photons with respect to the black 
body seed photons directly visible. The BMC model repro- 
duced rather well the data {x^—0.93, 64 d.o.f.) both at high 
and low energies where it was only modified by photoelec- 
tric absorption and a gaussian Une (see Fig. 4.). The best 
fit parameters are kT66=1.40±0.07 keV, a=1.37±0.9 and 
log(A)=-1.04l|!, g. We note that although the intercalibra- 
tion costant between XMM-Newton and IBIS is expected to 
be around 1, we found a smaller value. In principle this could 
imply some variability in the source flux between the obser- 
vations, however the most reasonable explanation is that the 
soft X-ray flux measured by XMM-Newton at apastron pas- 
sage represents the lowest quiescent state of the source while 
the hard X-ray flux, measured by IBIS over a much longer 
time interval, represents the average quiescent state along 
the eccentric orbit. 



4 ALTERNATIVE ESTIMATION OF THE 
NEUTRON STAR MAGNETIC FIELD 

As previously reported in section 2 and 3, IGR J18483— 0311 
has been detected in quiescence by XMM-Newton and IBIS 
with a luminosity in the range lO^^-lO^* erg s~^ and spectral 
shape similar to that during outburst activity. Such informa- 
tion could bo used within the framework of accretion from a 
spherically symmetric stellar wind to obtain an alternative 
and independent estimation of the magnetic field. In fact, 
the material in the background wind will be accreted onto 



the neutron star only if the magnetospheric radius _Rm is less 
than the corotation radius Rn (Rn > Rm), otherwise cen- 
trifugal forces will expel the matter (Davidson & Ostriker 
1973; Stella et al. 1983). We are aware that the above as- 
sumed model and propeller effect arc based on various and 
stronlgy simplified assumptions, however this is still a rea- 
sonable approach to obtain a rough and reliable estimate for 
the magnetic field strength and the eccentricity. The mag- 
netospheric radius is obtained by balancing the matter pres- 
sure to the magnetic field pressure, i.e. pV^ = B( R^^) /Sir 
where B{Rm) = Bns-RnsZ-^m (-Rns is the radius of the neu- 
tron star, Mns is the mass and Bns its surface magnetic 
field). Although the measured Mns in accreting X-ray sys- 
tems is known to be roughly in the range 1-2.5 Mq, neutron 
stars masses determinated at the highest accuracy display a 
normal distribution centered at the value of 1.35 M© with a 
very small dispersion of ±0.04 Mq (Thorsett & Chakrabarty 
1999). Bearing this in mind, here and in the following calcu- 
lations we assume the mass of the neutron star as a costant 
equal to its canonical value of 1.4 Mq. The corotation radius 
is given by: 



Rn = 



I 47r2 



1/3 



(1) 



where P is the spin period. Assuming a typical neutron 
star mass of 1.4 Mq and P=21.05 s, as measured for 
IGR J18483-0311, we get Rn=1.27xl0^ cm. 

The observed quiescent X-ray luminosities could be pro- 
duced only if the background wind material can reach the 
surface of the neutron star and impact onto it. However, the 
requirement of -Rn > Rm imposes important constraints. 

We calculated the Bondi-Hoyle X-ray luminosity at or- 
bital phases (p = and (p = 0.52 assuming a supergiant mass 
Mob = 33 M© and a radius Rob ~ 33.8 R© (Scarlc et al. 
2008), and Mns = 1.4 Mq, i?NS = 10 km for the neutron 
star. Then we assumed the typical wind parameters for a 
star with same spectral type of the IGR J18483— 0311 super- 
giant: a mass loss rate M in the range 0.4— 3x 10~^ M© yr^^, 
the wind velocity law v{r) — Vao{l — -Rob/'")'', where 
Voo = 1200 - 1900 km s"\ /3 = 0.8 - 1.6 and we leave the 
magnetic field (B) and the orbital eccentricity (e) as free 
parameters. We calculated the expected X-ray luminosity 
of the neutron star at periastron and at the orbital phase 
<f) ^ 0.52, assuming that the X-ray emission is due to accre- 
tion of matter from the spherically symmetric stellar wind 
emitted by the supergiant companion. The neutron star ac- 
cretes only the fraction of the wind ejected by the OB star 
within a distance from the neutron star smaller than: 



2GMns 



here 7?a is called accretion radius, and Vre\ is given by: 
Vrelir) = (Wwind(r) - Vr{r)f -\- w^(r) 



(2) 



(3) 



where Vrei is the relative velocity of the wind material Vwind 
with respect to the radial Vr and tangential V4, components of 
the neutron star orbital velocity, and Cs is the sound velocity. 
The fraction of stellar wind captured by the neutron star is 
given by the mass ffux which goes through the area TfR^, 
and is given by: 



Mace = Pwindir)Vrelir)TrR^ 



(4) 



Hint of a cyclotron emission feature in IGR J18483—0311 5 



where Pwind('") is the wind density related to the rate of 
mass loss by the stellar wind (M) by means of the continuity 
equation. From Equations (2), (3), (4), we computed the X- 
ray luminosity Lx of a neutron star: 

GMns ,> (GMns)' 47rpw(r) 

where G is the gravitational constant. The allowed range 
of parameters (e and B) which reproduce the X-ray lumi- 
nosity at the periastron {(j> ~ 0) and apastron (0 ~ 0.52), 
calculated according to Equation (5), are B^6xl0" G and 
0.39<e<0.62. We note that the above estimated magnetic 
field is in very good agreement with the value independently 
inferred by us from the putative electron cyclotron line at 
~3.3 keV (~3-4xl0" G). 



5 SUMMARY AND CONCLUSIONS 

In this work we studied the SFXT IGR J18483-0311 in qui- 
escence using archival INTEGRAL and XMM-Newton data. 
As results, we report for the first time the spectral proper- 
ties (r=2.5±0.8) and luminosity measurement (1.32x10'^* 
erg s~^) of the average quiescent hard X-ray emission above 
18 keV, accumulated over a long time interval along the ec- 
centric orbit. On the contrary, the 0.5-10 keV observation 
took place in coincidence with the apastron passage and 
allowed the measurements of the lowest quiescent state of 
the source (1.3x10^^ erg s~^); its spectrum is best fitted 
by an absorbed black body with parameters consistent with 
those previously reported by Romano et al. (2009) from a 
Swift /XRT monitoring campaign. 

One important point deserves further and deeper atten- 
tion. An intriguing spectral emission line has been observed 
at ~3.3 keV in the XMM-Newton quiescent spectrum of 
IGR J18483— 0311 at apastron passage and we speculated 
about its physical origin. An explanation in terms of atomic 
emission lines (i.e. highly ionized Argon) appears unlikely. 
We attempt to ascribe it to an electron cyclotron emission 
line from an X-ray pulsar accreting at a low rate, as pre- 
dicted by Nelson et al. (1995,1993). Specifically, these au- 
thors predicted the possible detection of electron cyclotron 
emission lines in the X-ray spectra of magnetized and tran- 
sient X-ray pulsars during their low luminosity quiescent 
state (i.e <10^'^ erg s"'^). The energy line center is ex- 
pected to peak at energies in the range ~2-20 keV and it 
should be superposed on the underlying soft thermal emis- 
sion. Emission-like features, similar to the ones predicted, 
have been observed in the X-ray spectrum of a handful of 
transient X-ray pulsars in HMXB systems during low lumi- 
nosity states (Nelson et al. 1995). We point out that the 
emission feature observed from IGR J18483— 0311 could be 
similar to the one predicted by Nelson et al. (1995, 1993). 
In our specific case, the lowest quiescent X-ray luminosity 
state achieved during the apastron passage as well as the 
sufficiently long XMM-Newton observation could have pos- 
sibly favoured the detection of a putative cyclotron emis- 
sion line. This interpretation would imply a neutron star 
magnetic field value in the range (3 4)xl0^^ G, which is in 
very good agreement with that independently inferred by 
us in the framework of accretion from a spherically sym- 
metric stellar wind model for the source. It is also worth 



pointing out that the estimated magnetic field strength is of 
the order expected for neutron stars known to display X-ray 
pulsations in high mass X-ray binaries, i.e. ~ (0.5-4) x lO'^^ 
G. Consecutively, it is strong enough to affect the accretion 
flow within a few neutron star radii and cause the observed 
X-ray pulsations in IGR J18483-0311. 

To date, IGR J18483-0311 is the only SFXT, among ~ 
10 firm members of this class, which exhibits such intriguing 
emission feature with important implications on theoretical 
models. However, this could be explained with the lack of 
sufficiently long X-ray observations of SFXTs at apastron 
passage (i.e. during their lowest quiescent X-ray states) with 
appropriate and sensitive enough X-ray facilities. 

Unfortunately, we can go no further on these issues be- 
cause the available exposure time and statistics of the data 
prevent us from a more detailed investigation. Longer X-ray 
observations of IGR J18483-0311 using XMM-Newton, for 
example, are strongly needed in order to achieve a higher 
signal to noise. This would allow us a much deeper investi- 
gation, in order to support or reject our proposed interpre- 
tation in terms of electron cyclotron emission line. 
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